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ABSTRACT 

The Polaris Flare cloud region contains a great deal of extended emission. It is at high declination and high Galactic latitude. It was previously 
seen strongly in IRAS Cirrus emission at 100 microns. We have detected it with both PACS and SPIRE on Herschel. We see filamentary and 
low-level structure. We identify the five densest cores within this structure. We present the results of a temperature, mass and density analysis of 
these cores. We compare their observed masses to their virial masses, and see that in all cases the observed masses lie close to the lower end of the 
range of estimated virial masses. Therefore, we cannot say whether they are gravitationally bound prestellar cores. Nevertheless, these are the best 
candidates to be potential prestellar cores in the Polaris cloud region. 
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1. Introduction 

In this paper we present observations, performed with the ESA 
Herschel Space Observatory (Pilbratt et al 2010), of the Polaris 
Flare region. In particular we use the large collecting area and 
powerful science payload of Herschel to perform imaging pho- 
tometry using the PACS (Poglitsch et al 2010) and SPIRE 
(Griffin et al 2010) instruments. These observations were carried 
out as part of the guaranteed-time key programme to map most 
of the Gould Belt star-forming regions with Herschel (Andre et 
al 2010). The Polaris Flare was first detected in HI as a spur 
of gas that appears to rise more than 30° out of the Galactic 
Plane. This region is an area rich in IRAS cirrus emission (e.g. 
Low et al 1984), and is sometimes known as the Polaris Cirrus 
Cloud. It was mapped in CO by Heithausen & Thaddeus (1990). 
On the large scale this cloud appears to merge with the Cepheus 
Flare cloud (e.g. Kirk et al 2009), and both clouds extend to high 
Galactic latitude. 

One of the denser regions in the cloud is known as molecu- 
lar cloud 123.5+24.9, or MCLD 123.5+24.9 (e.g. Bensch et al 
2003) - hereafter MCLD 123 - at a distance of 150 pc (Bensch 
et al 2003). It shows strong extended IRAS 100-yum emission 
and is generally believed to be gravitationally unbound with a 
mass of -18-32 M (Grossmann et al 1990; Bensch et al 2003). 
A CO study by Falgarone et al (1998) revealed a curved filament 
in MCLD 123 in 13 CO and C 18 - both in the J=2-l transition. 
This filament is also apparent in some narrow velocity channels 
in the same transition of 12 CO (Falgarone et al 1998). 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



There is one IRAS source in the region, IRAS 01432+8725. 
This is listed in the IRAS catalogue as having a flux density at 
100 pm of 2.88 Jy, but only upper limits at the other IRAS wave- 
bands. There is also one Spitzer source that was only detected at 
a wavelength of 24 pm at coordinates R.A. (2000) = 01 h 58 m 
27.5 s , Dec. (2000) = +87° 40' 07" . It has a peak flux density 
at 24 pm of 1.3 mJy/beam, where the Spitzer beam at this wave- 
length is 7 arcsec. This detection lies in a Spitzer calibration field 
in an unpublished archival dataset (AOR 33136386). 



2. Observations 

The SPIRE/PACS parallel-mode Science Demonstration 
Observations of the Polaris cloud were performed on 2009 
October 23 (Operation Day 162) at wavelengths of 70 pm and 
160 pm with PACS, and at 250 pm, 350 pm and 500 pm with 
SPIRE. The 70- and 160-//m ~6deg 2 scan map was taken with 
60 arcsec/sec scanning speed. The field was observed twice with 
both instruments by performing cross-linked scans in two nearly 
orthogonal scan directions. The combination of nominal and 
orthogonal coverages reduced the effects of 1/f noise and better 
preserved spatial resolution. The SPIRE data were reduced 
using HIPE version 2.0 and the pipeline scripts delivered with 
this version. These scripts were modified, e.g. observations that 
were taken during the turnaround of the satellite were included. 
A median baseline (HIPE default) was applied to the maps and 
the 'naive mapper' was used for map making. 

The PACS data were reduced with HIPE 3.0.455 provided by 
the Herschel Science Center (HSC). We used file version 1 flat- 
fielding and responsivity in the calibration tree, instead of the 
built-in version 3. Therefore the error in the final reduced flux 
scale was corrected manually with the corresponding correction 
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Fig. 1. The densest part of the Polaris Flare region at some of the observed wavebands. Upper row: 160yum from PACS, and 250/mi 
and 350/mi from SPIRE. Lower row: False-colour image (where 160yum is shown in blue, 250/mi is shown in green, and 350/mi is 
shown in red), column density map (where red is <4, blue is 4-8, and yellow is >8 x 10 21 cm" 2 ), and colour temperature map (where 
blue is 10-11 K and yellow is 12-13 K). The contour levels on the column density map start at 4 x 10 21 cm -2 , and the interval 
between successive contours is 1.5 x 10 21 cm -2 . The same contours are repeated on the temperature map for ease of location. Five 
sources are seen above a column density of 4 x 10 21 cm -2 . These are labelled cores 1-5 (in order of increasing R.A.) on the last two 
panels and are discussed in the text. The loop (loop 1) discussed in the text (containing cores 4 & 5) is clearly visible in all images. 
The reddest features on the false-colour image are the coldest, and the loop shows up clearly as redder than the surroundings. 
Likewise in the temperature map, the loop shows up as blue, indicating that it is the coldest feature on the map. The position of the 
IRAS source (IRAS 01432+8725) is marked with a star on the last two panels (adjacent to core 4). 



values in the PACS wavelengths. Standard steps of the default 
pipeline were applied for data reduction starting from (level 0) 
raw data. Multi-resolution median transform (MMT) deglitching 
and second order deglitching were also applied. Baselines were 
subtracted from the level 1 data by high-pass filtering with a -1° 
filter width, avoiding obvious sources, whilst the full leg length 
was 2.5° in the parallel mode. 



The PACS data of this field include transients of unknown 
origin after each calibration block, which seriously affected the 
ensuing frames. We processed these observations using data- 
masking and a narrower high-pass filter width than the image 
size in order to mitigate the calibration block artifacts. In this 
process, we may have removed spatial scales larger than the fil- 
ter widths. The final PACS maps were created using the HIPE 
'MADmap' mapping method projected to the 3.2 and 6.4 arc- 
sec/pixel size for 70 and 160 /mi data, respectively. 



3. Results 

The Polaris Flare dark cloud region was observed at five wave- 
lengths - 70, 160, 250, 350 and 500/zm. Figure 1 shows some 
of the main results. Only the densest part of the mapped region 
is shown. The area shown is just over half a degree square. The 
upper row of Fig. 1 shows the data from three of the wavebands: 
160/im from PACS; and 250/im and 350yum from SPIRE. The 
data have been smoothed to a common resolution of 24 arcsec, 
the approximate resolution of the 350-yum data. The images have 
also been re-gridded onto 10x10 arcsec pixels. 

The lower row of Fig. 1 shows some images derived from 
the raw data: a false-colour image; a column density map; and 
a colour temperature map. The contours on the column density 
map are at 4, 5.5, and 7 x 10 21 cm -2 . These are repeated on the 
temperature map to assist in source location. The Polaris cloud 
is clearly seen, and the raw data show a complex structure that 
is broadly similar at all wavebands. There are a number of fil- 
amentary structures seen in the data, with a few brighter cores 
embedded in the cloud. 
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Table 1. The physical properties of the cores. The Galactic latitude and longitude, as well as Right Ascension and Declination, are 
listed, along with the assumed distance. The radius of each core was measured by taking the column density map in Fig. 1 and 
measuring the equivalent radius of the contour that encircled a column density of 4 x 10 21 cm -2 , and this radius is given in pc. 
The full- width at half maximum (FWHM) is the geometric mean FWHM measured on the peak of each source. The integrated flux 
density within this contour at each of the Herschel wavelengths is listed in Jy. The absolute uncertainty in the flux densities is ±15%. 
As described in the text, a temperature, peak column density and mass were derived, and these are also listed. The uncertainty in the 
masses could be as high as a factor of 2. A mean volume density is given, assuming each core is spherical, within the given radius. 
Finally, a virial mass for each core is estimated using CO and HCO+ linewidths. 



Parameter 


Core 1 


Core 2 


Core 3 


Core 4 


Core 5 


Galactic longitude (2000) 


123.388 


123.511 


123.559 


123.687 


123.690 


Galactic latitude (2000) 


+24.928 


+24.915 


+24.856 


+24.894 


+24.931 


Right Ascension (2000) 


01 h 34 m 01.9 s 


01 h 44 m 51.6 s 


01 h 47 m 40.8 s 


01 h 59 m 42.7 s 


02 h 00 m 58.7 s 


Declination (2000) 


+87° 45' 42" 


+87° 43' 35" 


+87° 39' 33" 


+87° 39' 53" 


+87° 41' 58" 


Distance (pc) 


150 


150 


150 


150 


150 


Radius (pc) 


0.023 


0.035 


0.032 


0.035 


0.034 


FWHM (pc) 


0.023 


0.039 


0.027 


0.042 


0.038 


^cU(Jy) 


<0.18 


<0.18 


<0.18 


<0.18 


<0.18 


K^w 


4.26 + 0.07 


10.19 + 0.07 


7.56 + 0.07 


7.63 + 0.07 


5.44 + 0.07 


FS^(Jy) 


6.74 + 0.04 


16.77 + 0.04 


13.10 + 0.04 


15.35 + 0.04 


13.26 + 0.04 


FS^(Jy) 


3.61 + 0.02 


8.98 + 0.02 


7.27 + 0.02 


9.05 + 0.02 


8.50 + 0.02 


F ?oWW 


1.72 + 0.02 


4.33 + 0.02 


3.50 + 0.02 


4.53 + 0.02 


4.40 + 0.02 


Temperature (K) 


12+ 1 


12+ 1 


12+ 1 


11 + 1 


10 + 1 


N(H 2 )^(xl0 21 crrr 2 ) 


6 + 3 


7 + 3 


9 + 4 


13 + 5 


13 + 5 


Mass (Mq) 


0.1 


0.3 


0.3 


0.4 


0.5 


n(H 2 )mean (cm -3 ) 


~5xl0 4 


~4xl0 4 


~4xl0 4 


~5xl0 4 


~7xl0 4 


M vir (Mo) 


-0.3-0.5 


-1.0-1.5 


-1.0-1.5 


-1.0-1.5 


-1.0-1.5 



There is a filamentary loop seen in all images that is cen- 
tred roughly at Galactic coordinates 1 = 123.67, b = +24.89 - 
R.A.(2000) = 01 h 58 m , Dec.(2000) = 87° 40'. We here label this 
feature loop 1 . This is the same curved filament as was seen by 
Falgarone et al (1998) in 13 CO. They interpreted this as an edge 
of a cloud core. However, in the continuum we see it is clearly a 
loop with no filled centre. It was also detected in various transi- 
tions by Grossman & Heithausen (1992). 

There is also a filament with an apparent bifurcation 
at roughly Galactic coordinates 1 = 123.48, b = +24.90 - 
R.A.(2000) = 01 h 42 m , Dec.(2000) = 87° 43'. A bright core re- 
gion is seen at the head of this bifurcation, which may be bro- 
ken up into three components in the 160-//m data. The mean 
off-source pixel-by-pixel 1-cr variation on the N(H2) map varies 
from 1.2 to 1.5 x 10 21 cm" 2 . Hence, we adopt a value of 4 x 10 21 
cm -2 for the 3-cr contour. 

Five sources are seen in the column density map above a 
column density of 4 x 10 21 cm" 2 . We here label these cores 1-5 
in order of increasing Galactic longitude - see lower right panels 
of Fig. 1. We list the core positions and their assumed distances 
in Table 1 . The core mentioned above at the bifurcated filament 
is core 2, and loop 1 contains cores 4 and 5. 

The IRAS source IRAS 01432+8725 lies an arcminute to 
the west of core 4. We believe this offset is sufficient that the 
two sources are different (the IRAS FWHM at 100 /mi is 44 arc- 
sec). Therefore, none of the cores is associated with an infrared 
source, and so these are all candidate starless cores (Myers et 
al., 1987). The IRAS source is coincident with the centre of the 
loop, and may in fact be loop 1 itself, as IRAS point sources that 
only show up at 100 jim have often in the past been shown to be 
simply bits of cirrus. The Spitzer source may be foreground, as 
it is only seen at the shortest wavelengths. 

The reddest features on the false-colour image are the cold- 
est, and loop 1 shows up clearly as redder than the surroundings. 
Likewise in the temperature map, the loop shows up as blue, in- 



dicating that it is the coldest feature on the map. Cores 4 & 5 
appear to be the densest features on the map, with peak column 
densities in excess of 10 22 cm -2 . The column density contour of 
4 x 10 21 was selected as the core boundary in each case. The 
radial sizes of the cores were estimated from the images as the 
equivalent radius of a circle with an area equal to that contained 
by the core boundary. The derived equivalent radii are listed in 
Table 1 . Flux densities were measured within the core boundary 
contour in each case, and these are also listed in Table 1 . 



4. Core properties 

The core properties were estimated from the maps of column 
density and temperature. The flux densities of the pixels co- 
incident with the column density peaks of each core are plot- 
ted on the spectral energy distributions (SEDs) shown in Fig.O 
Modified blackbody curves were fitted to the flux densities, and 
these are also shown in Fig. [2l These are the same fits that 
were used, pixel-by-pixel, to construct the column density and 
temperature maps shown in Fig. [TJ The form of the fit (c.f. 
Hildebrand 1983) that was used in each case is 



:Q£ y (r)(m H yuN(H 2 K), 



(1) 



where F v is the flux density at frequency v, Q is the solid angle 
of each pixel, B V (T) is the blackbody function at temperature T, 
mnii is the mean particle mass (mu is the mass of a hydrogen 
atom and [i was taken to be 2.86, assuming the gas is -70% 
H2 by mass), N(H2) is the column number density of molecular 
hydrogen, and k v is the dust mass opacity. 

We used the pixel-by-pixel SED fits to calculate the col- 
umn densities, and hence the core masses. The value of k v that 
should be used has been the subject of much controversy. Here 
we adopt the dust opacity recommended by Henning et al (1995) 
and Preibisch et al (1993) for clouds of intermediate density - 
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Fig. 2. Spectral energy distributions of cores 1 to 5. The peak flux density in a single 10 x 10 arcsec pixel was measured. This is 
shown on a log-log plot of AS ^ versus A. The data are shown with 15% uncertainty error-bars. The upper limits at 70 yum are shown 
as arrows. The solid lines are grey-body fits of the form described in the text. The temperatures of the fits are listed in Table 1. 



n(H2) < 10 5 cm 3 - and we assume a standard gas to dust mass 5. Conclusions 



ratio of 100. This is a similar parameterization of the dust opac- 
ity to that used by Beckwith et al (1990), namely that 



k v = 0.1cm 2 g" 1 x (v/lOOOGhzy 5 , 



(2) 



where we have set the dust opacity index J3 to be equal to 2. 

This is also consistent with the value used by Andre et al 
(1993; 1996) and by Kirk et al (2005) for prestellar and star- 
less cores. The peak column densities and the temperature at the 
peak are listed in Table 1 . The mass of each core was calculated 
by integrating the column density map within the selected core 
boundary. This is also listed in Table 1 . From these, the volume 
densities were calculated, assuming that the cores are spherical. 
These, too, are listed in Table 1. 

The temperatures of the SEDs are listed in Table 1 . These are 
all quite low, with values of 10-12 K. This, and the lack of NIR 
emission from the cores, implies that the star formation process 
has yet to begin within these particular cores. Hence they are 
starless cores. This means that these cores should have no inter- 
nal heating and should be heated solely by the external radiation 
field. This is similar to what is seen in other low-mass starless 
and prestellar cores (e.g. Ward-Thompson et al 2002; or for a 
review see Ward-Thompson et al 2007). 

Core 5 was observed in the submillimetre by Bernard et al 
(1999). Our results are consistent with their findings, allowing 
for the very different resolutions of the two sets of observa- 
tions. Falgarone et al (2009) measured the mean CO linewidths 
in MCLD 123. For so-called 'bright' regions (i.e. high column 
densities) they found a mean linewidth of 0.4 kms" 1 in this re- 
gion. Heithausen et al. (2008) observed MCLD 123 in a number 
of transitions and found mean linewidths from 0.2 to 0.4 kms" 1 . 
Using this range of values we estimated a range of values for the 
virial masses of the five cores, and list these in Table 1. 

We note that all of the cores have masses that are below the 
virial masses that we have estimated. However, given the uncer- 
tainties in the mass calculations, they could be consistent with 
the lower limit of the range of virial masses. Hence, we can 
only say that they may or may not be gravitationally bound, and 
these 5 cores may be on the edge of possibly becoming prestellar 
cores. Note that this is very different from the cores found in the 
Aquila region (Andre et al, 2010), a large fraction of which are 
clearly gravitationally bound - c.f. Fig. 4 of Andre et al (2010). 
Nevertheless, bound or unbound, the 5 cores we have selected 
are the closest to being gravitationally bound of any of the star- 
less cores in Polaris. 



We have presented Herschel data of the Polaris Flare dark cloud 
region, and in particular the region MCLD 123. We found a great 
deal of extended emission at wavelengths from 70 to 500 jim 
with both PACS and SPIRE. We noted some filamentary and 
low-level structure. We identified the five densest cores within 
this structure. We carried out a temperature, mass and density 
analysis of the cores. We compared their observed masses to 
their virial masses, and found that the observed masses are on 
the lower limit of the range of their estimated virial masses, and 
thus we cannot say for certain whether they are gravitationally 
bound. 
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